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Abstract. The study of the spin structure of the nucleon and of the effects rising from the 
quarks transverse momentum are part of the scientific program of COMPASS, a fixed target 
experiment at the CERN SPS. The azimuthal asymmetries which appear in the cross-section 
of SIDIS off unpolarized targets can give insights on the intrinsic momentum structure of the 
nucleon and on the possible correlation between transverse spin and transverse momentum 
of the quarks. Here we present the new results for these asymmetries obtained from the 
COMPASS data collected with a 160 GeV/c positive muon beam impinging on a ^LiD target. 
The asymmetries are measured for both positive and negative hadrons, and their dependence 
on several kinematical variable has been studied. 



1. Introduction 

The SIDIS cross section for fongitudinally polarized leptons over unpolarized targets can be 
expressed [1] as 

do- (1 + (1 - yf) A / N Aim 

d.d,d.dj„.FjdPj = 2 ^o(l + ^.feMS'*.cos*.+ 

where Aq = Fuu and A[:^^^ = of reference [1]. The angle 4>h is the angle between the 

hadron production plane and the lepton plane, as described in figure 1, is the longitudinal 
polarization of the incident lepton and the quantities are given by: 
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The scaling variable is x = -p:^, where P is the target nucleon 4- vector, q = I — I is the 4 
momentum transfer (with I and /' the 4 momenta of the incoming and scattered lepton) and 
Q"^ = —q'^- y = ^ is the fraction of energy of the incoming nucleon transferred to the target 
nucleon. 

The sin (ph modulation has no clear interpretation in terms of Parton Model. Its amplitude 
is proportional to the beam polarization A; and it has been found to be small but different from 



Figure 1. SIDIS in the 
Gamma Nucleon System: 
P'* is the momentum of the 
produced hadron and 
its transverse component 
with respect to the virtual 
photon direction. 

zero (5% with a proton target) by the CLAS collaboration [2] which used a 6 GeV electron 
beam. 

The amplitudes of the cos (ph and the cos 2<j)h modulations have been measured in unpolarized 
SIDIS on a proton target by previous experiments and results have been published by EMC [3] 
and E665 [4] in a kinematical region similar to COMPASS and at higher energies by the ZEUS 
experiment [5]. More recent results, for positive and negative hadrons separately, have been 
presented by HERMES [6] and CLAS [13]. The origin of these two modulations is related with 
the so called Cahn effect, with the Boer-Mulders Transverse Momentum Dependent Parton 
Distribution Function (TMD PDF) and with the perturbative QCD corrections, which are 
expected to be important only for hadron transverse momentum Pji > 1 GeV/c. 

The contribution to the cos (ph and the cos Icph azimuthal modulation given by the Cahn effect 
comes from kinematics. It is calculated starting from the elastic quark-lepton cross section [7], 
taking into account the intrinsic transverse momentum of the quark inside the nucleon ky. 
The resulting cos <j)h amplitude is proportional to < > /Q and the cos 2(/)h amplitude is 
proportional to < > /Q^. Extraction of the mean value of k'^ from the measured amplitude 
of the cos (j)h modulation has been already presented [8] . 

The Boer-Mulders TMD PDF, which describes the correlation between transverse momentum 
and transverse polarization of a quark inside an unpolarized nucleon, has been introduced 
recently [9]. It contributes to the cos(j)h and the cos2(/)/j azimuthal modulation in unpolarized 
SIDIS where it appears convoluted with the Collins Fragmentation Function (FF). In some 
model calculations [10] it is predicted that the contribution of the Boer-Mulders function to 
the cos 2(j)h asymmetry is comparable with the Cahn effect and that it significantly differs for 
positive or negative hadrons. A first extraction of the Boer-Mulders PDF from previous data 
has already been preformed [11]. 

New results on these azimuthal asymmetries have been produced by the COMPASS 
experiment and they are shown here for the first time. 

2. Experiment 

The COMPASS experiment [12] is a fixed target experiment at SPS at CERN. It uses a 160 
GeV/c /i"*" beam, with a natural polarization of —80%, and a solid state polarized target. From 
2002 to 2006, data has been taken with a polarized ^LiD target. The targets were either 
transversely or longitudinally polarized to accomplish the wide experimental program on the 
spin structure of the nucleons. Up to 2004 the target consisted of two separate cells filled with 
the same material, oppositely polarized. The polarization of the cells have been reversed every 
week, during the data taking periods, in order to minimize systematic effects. 

3. Analysis 

The data used to extract the results are the part of the 2004 run in which the target was 
transversely polarized. To extract the amplitudes of the azimuthal modulations present in 
equation 1, the data with opposite polarization have been combined to cancel the small spin 




dependent contributions. Only the events with the /U^ scattered in the downstream cell have 
been used in order to have a larger angular acceptance for the detected hadrons. 

Apart for the almost standard DIS cuts: > 1 GeV^, 0.003 < x < 0.13, 0.1 < y < 0.9, 
W > 5GeV / (where W is the invariant mass of the 7*P system), a cut on the transverse 
hadron momentum has been introduced 0.1 < Pj^ < 1.0 GeV/c where the lower limit is chosen in 
order to have a good resolution on the azimuthal angle, while the upper limit ensures negligible 
pQCD corrections. A cut on z = EijE^* with 0.2 < z < 0.85 is applied to reject hadrons from 
the target fragmentation (lower limit) and away from possible diffractive effects at z = 1. Finally, 
cuts on the virtual photon angle (calculated with respect to the nominal beam direction), asking 
O^J^'* < 60 mrad, and y > 0.2 are applied for a better hadron acceptance. 

Dedicated Monte Carlo simulations have been performed to correct the hadron azimuthal 
distributions for the apparatus acceptance. For each of these simulations: the SIDIS event 
generation is performed with LEPTO; the interaction between particles and the experimental 
apparatus together with the detectors response is simulated with COMGEANT, a software based 
on GEANT 3; the Monte Carlo data are finally reconstructed with the same software (CORAL) 
used for the real data. 

Before evaluating the effects of the apparatus acceptance, extensive studies have been 
performed comparing the kinematical distributions of the reconstructed real and Monte Carlo 
data, in order to improve the agreement but limiting as much as possible the tuning of the 
LEPTO parameters. The ratios data/Monte Carlo for the adopted setting of the LEPTO 
parameters and after the standard DIS cuts are shown in figure 2 for the DIS variables and 
in figure 3 for the hadron variables. The agreement is good enough to correct the azimuthal 
distributions and to get the final results. 
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Figure 2. Ratios be- 
tween kinematical distribu- 
tions obtained from real 
data and Monte Carlo data 
for the DIS variables: x, y, 
Q2 and W. 



The azimuthal acceptance has been evaluated as: Acc{(ph) = ^rec{4''h^)/Ngen{4'h''^)^ where 
^rec{(l>h) is the reconstructed hadron distribution (after all analysis cuts) and Ngen{(ph) is the 
generated hadron distribution, obtained from Monte Carlo simulations. 

The hadron azimuthal distribution corrected for the apparatus acceptance Ncorr{4'h) = 
N{(j)h)/Acc{(j)h) has been fitted, binning alternatively in x, z and Pji, for positive and negative 
hadrons separately, with a four parameter function: Ncorr{4'h) = -^o • (1 -|- Pcos(</>h) " cos((/)/i) -|- 
Pcos(2<^h) ■ cos(20/i) -|- Psin(0h) " sin((/)/j)) where the free parameters pi are the amplitudes of the 
three independent azimuthal modulation (pcos^^ = ^i^m^^^)- 
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Figure 3. Ratios be- 
tween kinematical distribu- 
tions obtained from real 
data and Monte Carlo data 
for the hadron variables: 
the energy Ej^^i^ and the po- 
lar angle 9^^^, both calcu- 
lated in the Laboratory Sys- 
tem, and z. 



A typical hadron azimuthal distribution from real data N{(f)fi), the corresponding acceptance 
from the Monte Carlo simulation Acc{(j)h) and the corrected distribution Ncorri4'h) shown 
in figure 4, figure 5 and figure 6 respectively. 
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Figure 4. Measured azimuthal 
distribution N{(f)h)- 
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Figure 5. 
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Figure 6. Measured azimuthal distribu- 
tion corrected by the acceptance Ncorr- 
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The amplitudes of the azimuthal modulations introduced by the acceptance have been 
studied at length in the given kinematics. The azimuthal acceptances have been calculated 
in 2 dimensions where one variable is chosen among x, z or and the other among a set of 
other valuables on which the acceptance could have a relevant dependence like x, z, P^, y, 9l°i' 
and Q^. The kinematical regions in which the corrections are large have been identified and 
excluded; the most relevant cuts Ol"^' < 0.06rad and y > 0.2. The overall acceptance turned 
out to be flat. As an example the amplitude a^os rpiS'-"^ : Pt) the coscph modulation extracted 



from a fit to Acc{<^h) = Nrec{(l>h) / ^gen{(f>h) Calculated in each 2 dimensional bin in x and Pji is 
shown in figure 7. Similar pictures have been obtained studying other azimuthal modulations 



Figure 7. The amplitude 
dcos 4>h{^ J Pt) of the cos (f)h mod- 
ulation extracted from the fit 

to Ai(l)h) = NreciM/NgeniM 

calculated in each 2 dimensional 
bin in x and Pji- The absolute 
value of the amplitude is less than 
15% over the whole kinematical 
range and it is does not show any 
partcular trend. 

and for the combinations of the other kinematical variables. Thus the method used to extract 
the asymmetries should give the same results both using different event generators and with 
acceptances extracted in multi dimensional bins. 

4. Systematic studies 

The acceptance has been also calculated using two other Monte Carlo sets of data produced 
by changing the PDFs set and the LEPTO parameters relative to the transverse momentum 
distribution in the fragmentation process and the Lund Fragmentation Function. The agreement 
between the real data and these two sets of Monte Carlo data is worse than that of figure 2 
and figure 3, and shows opposite trends in the x, y and ^^J^f, distributions. For this reason they 
have been used to check the whole procedure. As an example of this test, in figure 8, A^!^^^ 
from positive hadrons extracted using the 3 different Monte Carlo is shown. As can be seen 
the differences are small as expected if the azimuthal acceptance is really flat in the kinematical 
variables. 

The difference of the asymmetries extracted using the 3 Monte Carlo has been used to evaluate 
the systematic error. To the systematic error contributes also the difference between the results 
obtained from the data taken with the transversely polarized target and the one obtained from 
the data taken with the longitudinally polarized target. In this second case the configuration of 
the apparatus is different, some detectors have been moved and many parameters changed and 
so, for the longitudinal target polarization data, a different Monte Carlo sample has been used to 
extract the unpolarized asymmetries. The effects due to possible inefficiencies of some detector 
planes have been also evaluated and found to be smaller than the previous ones. The sum in 
quadrature of all these three different contributions has been used to calculate the systematic 
error, which is twice the statistical one. 

Other tests have been performed finding negligible systematic effects. In particular the 
radiative correction has been calculated for in the given kinematics using RADGEN and are 
relatively small. The amplitudes of non physical modulations have been added in the fit and 
turned out to be compatible with zero. 

5. Results 

The results for the azimuthal asymmetries extracted from the COMPASS ^LiD data are shown 
in figure 10 for the A^^, amplitude, in figure 11 for the A^^2d> amplitude. 
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Figure 8. Compari- 
son between the A^^^^ 
azimuthal asymmetries 
extracted using 3 Monte 
Carlo samples using differ- 
ent event generations. The 
blue points are calculated 
with the Monte Carlo 
simulation used to extract 
the final results. 
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and in figure 12 for the A^^^^ amplitude. They are shown separately for positive (top) and 
negative (bottom) hadrons and are given as functions of x, z or PJp. 

The largest signal is given by the A^^^^ amphtude: up to 15% and with a strong dependence 
both on z and PJp. The dependence on x is less evident and its interpretation should take into 
account the correlation between x and Q^, shown in figure 9. The Al!^> asymmetry should be 

COS (pfi 

mainly due to the Cahn effect and the unexpected dependence of the results from the charge of 
the produced hadron could hint both to a different value of the intrinsic transverse momentum 
kT for the u and the d quarks or to the Boer-Mulders TMD PDF. This function appears here 
as an higher twist effect while comes at leading twist in the ^^^20^ &™plitude, thus is expected 
to give a relevant contribution to this last asymmetry. 

The A^^2(Pf^ asymmetry is smaller showing a strong trend as function of all the three studied 
kinematical variables x, z or P^. The dependence of A^^2(p^ hadron charge is clearly 

seen and could be related to the Boer-Mulders function and thus to an important role of the 
correlation between the intrinsic transverse momentum and spin inside the nucleon. 

A clear hadron charge dependence has been measured also by the HERMES Collaboration [6] 
both for the A^^^^ and the A^^^ij,^ amplitudes. Any comparison with the COMPASS results 
should take into account also for the different kinematical range covered by the two experiments. 

The azimuthal asymmetries show a positive signal for the positive hadrons while are 

compatible with zero for the negative ones. A signal of the same sign has been measured by the 
CLAS Collaboration [13] using a 6 GeV electron beam on a proton target. 

Summarizing, the COMPASS SIDIS measurements with unpolarized targets show large 
asymmetries and effects sometimes unexpected. Comparison with model calculations and the 
outcome of the other experiments will allow to clarify the present description of the nucleon 
structure including the role of the intrinsic transverse momentum and its correlation with the 
transverse spin. 
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Figure 10. Re- 
sults for the A^^^^^ 
azimuthal asymme- 
tries as functions 
of X, z and Py, 
obtained for posi- 
tive hadrons (up- 
per plot) and nega- 
tive hadrons (lower 
plot). Grey bands 
are the systematic 
error 
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Figure 11. 

Results for the 
^co^24>u azimuthal 
asymmetries as 
functions of x, 
z and P^, ob- 
tained for positive 
hadrons (upper 
plot) and negative 
hadrons (lower 
plot). Grey bands 
arc the systematic 
error 

[psys = 2 • (Jstat)- 
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Figure 12. Re- 
sults for the A^P^^ 
azimuthal asymme- 
tries as functions 
of X, z and Py, 
obtained for posi- 
tive hadrons (up- 
per plot) and nega- 
tive hadrons (lower 
plot). Grey bands 
are the systematic 
error 
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